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We perform various resummations of the hot QCD pressure based on the actual knowledge of the 
perturbation series which includes the ~ g^\n{l/gs) and part of the ~ gf terms. Resummations are 
performed separately for the short- and long-distance parts. The ~ pf term of the short-distance 
pressure is estimated on the basis on the known UV cutoff dependence of the long-distance part. 
The resummations are of the Pade and Borel-Pade type, using in addition the (Pade-)resummed 
expression for the squared Debye screening mass m| and, in some cases, even for the EQCD coupling 
parameter g^. The resummed results depend weakly on the yet unknown ~ gf terms and on the 
the short-range renormalization scale, at all temperatures. The dependence on the long-range 
^ , renormalization scale is appreciable at low temperatures T < 1 GeV. The resulting dependence of 

O ' pressure on temperature T is compatible with the results of the lattice calculations at low T. 
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I. INTRODUCTION AND DESCRIPTION OF THE APPROACH 



> 

(N 

^^ \ Today we have at our disposal a well-elaborated technique for perturbatively treating field theory at finite tem- 

perature and/or chemical potential 1]. This formalism goes far beyond the ordinary (T = 0) perturbation theory 
insofar as a correct treatment of infrared divergences (specifically those which are connected with T 7^ 0) requires 
f^ . partial resummation of infinitely many specific diagrams. This implies - among other things - that the infrared 
■^ ■ convergent perturbative expressions come out as series in g^ (the fundamental coupling constant) rather than in gl 

II ■ During the last decade, several physical quantities, the most prominent example being the free energy density F of 

Q^! the quark-gluon plasma, have been calculated within this formahsm up to 0{gl) |1,|31 and (partially) even to 0{g^) 
Ih ' SB Disappointingly, in spite of the relatively high orders available, the results are of very limited applicability even at 
very high temperatures (T > 1 TeV) where gs certainly becomes small. In fact, if successive terms in the perturbative 
series are added, the corresponding truncated sum changes dramatically, jumping up and down. Furthermore, the 
unphysical dependence on the renormalization scale is strong and seems to become even stronger with increasing 
order (the renormalization scale in MS scheme will be denoted as Jl, and in a general scheme as /z). Both effects 
C^ ' considerably reduce the reliability of the perturbative results for representing physical quantities and have therefore 
been object of intensive theoretical studies. 

Ways out of the convergence dilemma are in general looked for by either performing some clever resummations 
or by reorganizing perturbation theory in some way (actually, in some cases, this amounts to the same thing) [l|. 
Several specific approaches have been applied up to now. Among them are quasi-particle models [a|, ^-derivable 
approximations [3 , screened [3 or optimized |3] perturbation theory, hard thermal loop perturbation theory [l^ . 

Within the present paper we concentrate mainly on the question of renormalization scale (RS /I) dependence of the 
resummed results, in MS scheme, as well as on the dependence of the results on the yet unknown part of the gf-term in 
truncated perturbation series. The resummation approaches mentioned above do have some residual renormalization 
scale dependence. Furthermore, the conventional choice for the RS JI at finite T (/I ~ 2'kT) is not natural since 
different energy scales get involved in all calculations (see later). 

In order to improve this situation, we replace the (partially resummed) perturbation series by approximants which 
are more stable under the variation of Jl. The Fade approximants (PA's) [ill . Il2 . Ilil [lj| and other Pade-related 
resummations, such as Borel-Fade (BFA's) [lj| and modified Baker-Gammel approximants (mBGA's) jl6, 17], are 
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known in general to reduce the unphysical /I-dependence significantly (mBGA's even entirely). Usually these approx- 
imants are applied to a given truncated perturbation series (TPS) of a given order, and they fulfill in addition the 
"minimal" requirement: upon re-expanding them in powers of the coupling parameter, they reproduce the TPS to 
the given order. In this sense they are sometimes considered to be some sort of mathematical artifice without any 
deeper physical motivation. This opinion is delusive, however, because of several reasons: 

1. The basis of all field theoretic approaches is the path integral expression for the generating functional which has 
to be attributed a specific mathematical meaning. In ordinary perturbation theory this is achieved by Taylor 
expanding the corresponding exponential of the interaction Lagrangian leading to the well-known power series. 
But this is by no means better or more natural than by approximating it in terms of other approximants (e.g.. 
Fade) - we are only less used to it. In some sense Fade or Fade-related sequences seem even more adapted than 
power series since they allow for (pole) singularities, and we know that most physical amplitudes do include 
singularities. 

2. The asymptotically divergent nature of the expansion in powers of the coupling parameters represents a practical 
problem concerning the actual evaluation, and calls for, among other things, nonperturbative information in 
order to fix the renormalon ambiguities and to produce unique predictions. On the other hand, due to the 
generally better convergence of Fade and Fade- related approximants, one hopes that the convergence of the 
corresponding sequence of the (Fadc-)resummcd TFS's will improve and that it would converge to a specific 
prediction. The Fade sequence, in fact, with increasing order shows convergence under rather general conditions 
even when the corresponding TFS is (asymptotically) divergent |ll| . 

3. It can be shown that FA's (and mBGA's) can be represented by weighted averages of running coupling parameter 
(on e-loop running in the case of FA, exact running in the case of niBGA) at specific values of reference momenta 
[iM US > and are thus a realization of the fact that each physical quantity in field theory is characterized by 
specific values of the momentum flow. In this sense they represent different approximations to the Neubert's 
formalism [l9J of renormalization improved perturbation theory. 

From the diagrammatic point of view, FA and Fade-related methods represent a resummation of those (infinitely 
many) diagrams which - when added - approximately (or exactly) cancel the /I-dependence of the given quantity. 

FA's, BFA's and mBGA's have been applied to TFS for several QCD or QED quantities at zero temperature (see, 
for example, Refs. Jl3 . [l5l[l7| ) and have often led to significant improvement of the RS-dependence problem. Recently, 
same types of approximation have also been used for similar purposes in finite temperature gauge theory: FA's in 
Refs. 20, 21]; somewhat related Borel methods in Refs. 22], by discerning some information on renormalons (for 
renormalon properties in (p'^ theory at finite T, see Ref. |23]). The first FA resummations at finite-T [23] consisted in 
simply replacing the available power series (in powers of ^s) for the entire QCD free energy by various FA's. Although 
the results demonstrate a weakened T^-dependence, we regard some aspects of this approach as problematic. Our 
reservation is due to the simple observation that two ingredients - (a) the infrared stable TFS of thermal field theory, 
and (b) the Fade(-related) approximations applied to this TFS - imply resummation of infinitely many terms. The 
corresponding classes of diagrams are, however, neither equivalent nor disjunct. Therefore, care has to be taken to 
avoid double counting and to disentangle the various resummations. 

Recently, we have developed a formalism to consistently treat this problem pl|. The procedure is the following: 
Consider some physical quantity and its TFS (to a given order) in thermal field theory. We specifically have in mind 
the static pressure pqcb of the QCD plasma, i.e., the negative of the free energy density -Fqcd- It is connected to 
the partition function Z by the relation 

T 
PQCD = -Fqcd = ,lim TrlnZ, (1) 

nl/T 

■ j3 



Z= VA^^V^V^je^pl- dr d-'xCQCD] , (2) 

where £qcd is the (Euclidean) QCD-Lagrangian density, T is the temperature, V is the three-dimensional volume, 
and the renormalization convention PQcuiT = 0) = is taken; Qs is the QCD gauge coupling parameter. The 
corresponding TFS has the generic structure 

P - Pidcai [Co + C25? + C35f + C^gt + C59! + Cegt + ■■■]. (3) 

where the coefficients Ci may include contributions of order In^s. The terms proportional to odd powers of gg 
[or, equivalently, fractional powers of a = g'^ / {Att'^)] are exclusively due to resummation of specific (ring, or daisy) 



diagrams, but the latter also contribute to terms proportional to even powers of gs- Clearly gs and the coefficients 
Ci depend on the chosen renormalization scheme and, in particular, on the renormalization scale Jl (ex T). Due to 
the asymptotic freedom we expect that g^ becomes sufficiently small at large temperature T. In the following, we 
concentrate on such a large-T situation. 

The first step needed for the application of the formalism of Ref. [21| is the separation in Eq. ^ of the purely 
perturbative part from the contributions stemming from the (ring-)resummations. This can be done in an unambiguous 
and consistent way, since the resummation terms represent exclusively the contributions of the (bosonic) zero modes 
to the considered physical quantity. Note in this respect that resummation is needed purely for taming the finite-T 
infrared divergences, and only the infrared parts of the higher order diagrams contribute to the (finite) order terms. 
Therefore, the identification of the contributions from the zero modes (long-range contributions) is needed. This can be 
achieved either by integrating out directly the high momentum regime or, more elegantly, by utilizing the effective field 
formalism |3,|23 based on dimensional reduction method Refs. 3, 27, 28, 29]. The idea behind is that (at high enough 
temperature, with gs sufficiently small) the four-dimensional thermal QCD at length scales ^ \/{gsT) is equivalent to 
an effective three-dimensional field theory of (static) boson fields. This effective theory represents the physics of the 
zero bosonic modes and thus reproduces all static correlations of the original QCD at the aforementioned distances. 
Consequently, the contribution of the zero modes to the partition functions or to the pressure can be calculated 
by means of the effective theory. Within QCD the long-distance part can be further subdivided corresponding to 
the two regimes R ^ \/{gsT) (determined by color-electric screening) and R ~ l/{glT) (color-magnetic screening). 
Denoting the corresponding contributions to p (or F) by the subscripts M and G, respectively, the final decomposition 
corresponding to the three energy regions is 

P = P^+PM+PG ■ (4) 

Here, p-^ denotes the short-distance, and pM -I- pQ = pm+g the zero-mode long-distance contribution. 

Now comes one of the main points of the formalism of Ref. |2j| : Having the separation @ of the pressure into pe 
and pm+Gi we can argue that both are separately physical quantities, and thus both, when calculated to all orders, 
are separately /I-independent. The reason for it is shown by the following indirect argument: Consider any (static) 
two-point-correlator within the given field theory. It is a measurable quantity and so is its long-range behaviour. But 
the latter is determined solely by the exchange of zero modes (note that the exchange of a mode with frequency w„ 
contributes a term to the correlator which vanishes like exp(— |iu;„|i?) at i? — > oo, thus the only surviving contribution 
at large distances is the one with loq = 0). Thus the contribution of wq to each correlator is measurable (and therefore 
/i-independent) once we have fixed, by convention, the minimal value i?min of what we consider large distances 
(or: factorization scale Ae ^ 1/^min)- Since these correlators are derived from the partition function Z by means 
of functional derivatives, even the zero- mode contribution to Z, and thus also to p, is /I-independent (physical). 
Consequently, the remaining part of the measurable pressure - the short-range contribution pE represented by the 
ordinary perturbative terms - must also be /I-independent. The RS-independence of both parts separately has been 
demonstrated analytically in Ref. [2 if by using the known perturbation series (up to order g^). 

The separation and the /I-independence of the separate terms allows a more consistent application of Pade and 
Pade-related approximations: they are applied to TPS's of each quantity pE and pm-i-g separately."'^ We can illustrate 
the importance of using separate approximations with a simple example - the Pade-resummation of a quantity 
S = (Si -1-52), where 5i and 5*2 are separately physical quantities, each of them available as power series of a up to 
next-to-leading order (NLO): 

5j=a(l + r(,)a)+0(a3) (j = l,2). (5) 

When we apply to the TPS of the sum S a PA, say [1/1] (a), and then expand this back in powers of the coupling a, 
we obtain 



511/11 = 2 a 



1 



1- o(^(l)+''(2)) 



= 2 a 



1 + 7;{'i-(i)+r{2))a + -(r^i) + r^2) + 2r-(i)r(2))a2 



+ 0(a^). (6) 



The coefficient at ~ a^ here has a term 2r(i)r(2). Therefore, the Pade-resummed S contains mixing effects at ^ a?, 
i.e., an interference effect between the two amplitudes for the physical observables 5i and ^2. This is not acceptable, 
because S is the (incoherent) sum of S\ and 52- Therefore, the PA should not be applied to the entire sum 5, but 



^ Since only little information is available about pQ, we will find it more convenient to apply Pade (-related) resummations to Pm+G ^nd 
not separately to pM and PG- 



separately to Si and 5*2. This argument holds also when different PA's or Pade-related resummations are applied, 
and when the order of the TPS is higher. 

One of the consequences of the separate treatment of pE and pm+g is that the natural renormalization scale /I used 
in these two quantities should be of the order of the energy of the modes contributing to them: JI ^ 2ttT in pE, and 
71 < gsT inpM+G- 

Within the present paper we will apply Pade (PA) and Borel-Pade approximants (BPA's) to the evaluation of pqco- 
The modified Baker-Gammel approximants (mBGA's), which are /I-independent J^j or even renormalization scheme 
independent |l7|. are technically more involved. An analysis with mBGA's, which at least could serve as some kind 
of quality control of the general procedure, will be presented elsewhere. On the other hand, BPA's are also applied 
here and they represent an extension of the Pade analysis: PA's are applied to the corresponding TPS's of the Borcl 
transforms, and then the resummed quantity is obtained by Borel integration. This procedure gets its motivation 
from the hope that the Borel summation might defuse the notorious divergence problem of perturbation theory, as 
well as from the fact that the power expansion of the Borel transform has significantly better convergence properties 
than the original series and is thus more amenable to the Pade-type resummations. 

One technical remark should be added here - the separation of the energy range into a high and a low energy 
region requires introduction of a factorization scale A^; which defines a boundary between the two: 27rT > A^; > gsT. 
Consequently, the two contributions separated in the described way acquire an artificial dependence on Ae, although 
the sum of the two terms pe + Pm+g has to be AE-independent. Similarly, for the individual terms pe and pm+g to 
have physical meaning themselves, the factorization scale should be suitably chosen. Technically this means that the 
approximants for pE and pm+g, although independent of each other, should be chosen such that the A^j-dependence 
is minimized. 

In our previous paper |2l] | we applied this procedure to calculation of the free energy (pressure) both in QCD with nj 
(massless) quarks and in a ^''-theory. At that time the relevant TPS's had been calculated only up to terms of order 
g^ with the following implications for the perturbative structure of the relevant contributions: The short-distance 
term Fe {— Pe) is determined perturbatively up to NLO [in a = a{'p) = g^{jl) /{2Tr)'^] 

-pE/^ideal = PE/Pidcal = 1 -" B{nf)FE , (7) 

Fe = a{l + CE(n/,AE,7l)a} , (8) 



and so is the electric (Debye) screening mass ms 

~2 _ 1 1 2 

The long-distance part was a power series in gs [= gs (/!)] up to order g^ 



= a{l + Cuinf,^J.)a} (9) 



2 _ 

Fm = -PM = --^TrrV^pM , (10) 

OTT 

Pm = 1 + C'Mi(n/,AE,7l)gs + C'm2(«/)5s • (H) 

We refer to Ref. [211 for compilation of explicit expressions for B, Ce, Cm, Cmi, Cm2- These are rather short power 
series and they allow construction of only low order Pade approximants, e.g., PA [l/l](a) for Fe and TOe, and PA 
[0/2] (gs) for PM- In addition, the TPS for pu {— —Fm) is very strongly divergent, and the TPS's for Fe and for rh'^ 
are divergent to a somewhat lesser extent. Therefore, resummation results based on these TPS's should be taken with 
care. Nevertheless, the application of these low order approximants yielded results which were fairly stable under the 
variation of JI, although ~ at low temperatures T ^ 10 GeV - they deviated substantially from the lattice results 

[nmEi. 

In such a situation, an additional order in the perturbation series (i.e., terms of order g^) is much more than an 
additional tiny correction, but constitutes a significant enlargement of the basis for Pade approximations, since it 
adds additionally both to Fe and Fm and (for the first time non-vanishingly) to Fq. Therefore, the calculation of 
the ~ g^ \n{l/gs) contribution to the long-range part of the pressure published in Ref. jj| is extremely gratifying. 
The full 0{g^) contribution cannot be achieved perturbatively because of the well-known breakdown of perturbation 
theory due to incurable infrared divergences occurring at this order [Sfl, |3ll . What can be evaluated is the coefficient 
of the logarithmic ultraviolet divergence contained in Fq , because of the superrenormalizability of the corresponding 
effective three-dimensional field theory; and this is exactly the coefficient of the O (g^ ln(l/(7s))-term in Fq. The 
term purely proportional to g^ remains unspecified and has to be treated as a free parameter unless information from 
nonperturbative methods (e.g., lattice calculations) is inferred. 

Within the present paper we will utilize the new results |j| on the perturbation expansion of the QCD-pressure 
as much as possible in order to find approximants which are reasonably stable under RS-variation. This gives us 



the freedom of using higher order Pade or Borel-Pade approximants and thus treating the long-range part pm+g 
(= —Fm+g) in a more rehable way. Unfortunately, the full four-loop contribution to the short-range part p^ {— — Fg), 
which is in principle perfectly calculable within ordinary perturbation theory and would yield 0{a'^) correction to 
Eq. JSJ, is not yet available at the moment. Therefore, we have no direct basis for improving the approximants to p^. 
Nevertheless, we can obtain some restricted information about the 0{g^)-teims in pE from the requirement that pE 
be /I-independent and (pE +Pm+g) be AE-independent. A constant (jl- and AE-independent) term in the coefficient 
at 0{g^) in pE still remains unspecified, but its value can roughly be estimated. The new results of Ref. pl also allow 
a better approximation to the parameter g^ of the effective theory, which further improves the resulting predictions. 
In Sec.|n]we describe the separation of the QCD-pressure into contributions stemming from different energy regions 
and specify the known corresponding effective Lagrangians. We then present the available perturbation expansion of 
the long- and the short-range contributions and work out the effects of the factorization scales A^ and Am (which have 
not been explicitly disentangled in Ref. jj| because of their simplified treatment of the RS-dependence). In Sec. IIIII 
the short-distance and the long-distance parts of the pressure are resummed separately by Pade and/or Borel-Pade 
approximants, and the selection of approximants is narrowed down by requiring weak residual /I-dependence of pE and 
of pm+G: and weak A^j-dependence of (pE +Pm+g)- In Sec. lIVI the resummed results as a function of temperature T 
are presented, as well as arguments for narrowing down further the selection of acceptable approximants. In Sec. Ivl 
results of TPS evaluations as a function of T are presented, for comparison. In Sec. IVII our results are compared 
with the predictions obtained within other approaches, in particular with lattice results (which are available only for 
rather low temperatures), and finish with some concluding remarks. A short compilation of the Pade and Borel-Pade 
approximants is given in the Appendix. 

II. PERTURBATION EXPANSION OF LONG- AND SHORT-DISTANCE PRESSURE 

The basic information about the physics of a quark-giuon system in thermal equilibrium at temperature T is 
contained in the expression for the pressure pqcd(T) or, equivalently, the free energy density Fqcd(T). In d -I- 1 
dimensions (d = 3 — 2e) this is given by Eqs. CO-©; with (Px replaced by d'^x. Boundary conditions over the finite 
time T direction are periodic for bosons and anti-periodic for fermions. 

When the temperature T is above the masses niq of active quarks, and the QCD gauge coupling gg is small enough, 
there are three physically different scales involved: ~ 2ttT, ~ gsT, ^ g^T. The last two correspond to the color-electric 
and color-magnetic screening, respectively. The decomposition Q pqcd = Pe + Pm + Pg reflects the contributions 
from the aforementioned three energy regimes: pe are contributions from the modes with energies in the interval 
[AejOo], Pm from those in [Am, Ae], pq from those in [0,Am], where the factorization scales Ae and Am define the 
borders between the three energy regimes 

g^T < Am < gsT < Ae < 27rT . (12) 

The long-distance part pm +Pg = Pm+g is due to the bosonic zero (Matsubara) frequency mode, whereas pe contains 
all higher modes. For calculating the different parts analytically one most conveniently uses the method of effective 
Lagrangians ^, ^, .32J based on the dimensional reduction method ^ 23^ 2^ 2^ ■ Whereas the high-energy regime 
behavior, which is responsible for pE, is determined by the original {d + l)-dimensional QCD Lagrangian, the low- 
energy regime behavior, responsible for pm+Gj can be represented by a d-dimensional effective bosonic theory called 
electrostatic QCD (EQCD), such that 

PQCD = PE + ^lnfvAfVA^ex^p(~ J d'^xCEQCDJ , (13) 

/:eqcd = ^Tri^2 ^ r^^[^^^ ^^]2 ^ ^|TrA2 ^ ^W (TrAg)' + AJ^'^Tr^ + . . . , (14) 

where i = 1,.. .,d; Fij = {i/gE)[Di, Dj], Di = di — igE^i. The notation A^ = A'tT"' is used, with T" being the 
Hermitean generators of SU(3) normahzed to Tr T'^T''' = (l/2)(5"''. 

The four parameters of the effective theory are the Debye screening mass ttie (^ QsT), the coupling parameters g^ 
(~ g^T) and A*^^\ A^^-' (~ g^T) (the latter two are not independent if d = 3). These and the hard scale pressure pe 
(~ T^) are obtained by the well-known matching procedure: the parameters as functions oi gs,T (and the UV-cutoff 
Ae of vCeqcd) must be tuned in such a way that the effective theory reproduces the physical effects of the full theory 
at large distances. The dots in Eq. (|14|l stand for higher dimensional effective interaction terms which yield corrections 
6p ~ gjr^ and are thus neglected here. 

Due to the separate screening of the color-electric (mE ~ gsT) and the color-magnetic {ruM ~ 9^T) degrees 
of freedom within QCD, the above EQCD-action can be further decomposed into two parts corresponding to two 



physically different energy scales. The color-electric scales {'--^ QsT) are separated by integrating out A^ 

^ In / VAIVA-'o exp (- j d^'xCEQCuj = Pm{T) + ^ In / ^^" exp f~ f d^'xCuQCBJ (15) 

= Pm{T)+pg{T), (16) 



where 

-CMQCD^^Tr^^... (17) 

represents the Lagrangian density for magnetostatic QCD - the effective theory for the energy region below Am, 
typically energies ^ g'^T. Here Fij = (ilgy\)\Di, D.j\ with Di — di — i^M^i- 

There are two matching coefficients now: the long-distance [^ iKg^T)] pressure pM (~ ra^T ^ ggT^) and the 
coupling parameter g^ which is close to g^ [sF 



gl, = 5l (1 + 0{gl/mE)) = g| (1 + Oig, 



)) ■ (18) 



Now we will present the perturbation expansions for px (X — E, M, G) on the basis of the results of Ref. ;4i. The 
authors of Ref. 4] present each px/{T]l^'^'^) as^ an expansion in powers of the coupling parameters of the respective 
effective theory. We note that they used dimensional regularization in the MS scheme, thereby invoking a common 
renormalization scale JI for all px and for other matching coefficients. Therefore, their results do not contain explicitly 
the (physical) factorization scales Ag and Am, and contain the unphysical infinite terms ex 1/e which cancel in the 
sum pqcd = Pe-i-m+g- We will obtain our basic formulae on the basis of their aforementioned expansions by replacing 
in the logarithms of expansion coefficients their renormalization scale JI by the IR and/or UV cutoffs of the respective 
effective theory,^ and by removing all the terms proportional to 1/e in each px/{T'jj,~^^) and then taking e = in each 
Px- While the aforementioned 1/e terms cancel in the sum Pe+m+g / {TJi.^'^'^) , the effect of the common factor l//i~^^ 
also disappears (i.e., reduces to one) in this sum when e —^ 0. Which of the In/I-terms in the coefficients get replaced 
by logarithms of the UV scale and which by those of the IR scale of the effective theory - this is unambiguously 
determined by the requirement that the entire pe+m+g be independent of the factorization scales. On the other hand, 
the aforementioned procedure to eliminate first the 1/e infinities from each px/{T'jl~''^'^) and then replace e i— > to 
obtain px/T is certainly not unique. Other procedures would result in expansions for px which differ from ours by 
certain constant numbers (independent of scales) in some of the coefficients of the expansion in powers of the coupling 
parameters of the respective effective theory, this ambiguity being a manifestation of renormalization freedom. The 
overall sum pe-i-m-i-g would remain the same. We note that our procedure leads to the decomposition of Ref. Q, at 
least to the order available in that Reference. Further, we checked that the re-expansion in powers of QCD coupling 
parameter gs{p) of each obtained px reproduces for the sum pe+m+g the same expansion as the one obtained in 
Ref. 0. 



The expansion of pc starts at ^ g^ [^ g^) 

8-33 



PG = T56 



^(4^)4 



'saGlnfJ^] 



\2mM/ 



(19) 



43 157 o 

ac = TT^ w 0.195715 , (20) 

96 6144 ' ^ ' 

where the value l|^ was obtained in Q and tom = 'igl/i i'^ 5s ^) • The coefficient Sg denotes a dimensionless number 
which is not perturbatively calculable but is expected to be \Sg\ ^ 1 (note that Sac ~ 1-6). We will allow the 
following, rather generous, variation of Sq- 

-5 < Sg < +5 . (21) 

We note that, according to the procedure mentioned before, the result H19|l is obtained from the corresponding result 
for pg{T) / {T'jl~'^^) of Ref. 3 by removing the term 1/e, replacing in the logarithm Ji by the UV cutoff Am of MQCD, 
and then taking e — > 0. 



^ /i = /i(e^E/47r)i/2^ where fi is the renormaUzation scale in the MS scheme. 

3 For PG (MQCD), the UV cutofT is Am; for pM (EQCD), the IR cutoff is Am, and the UV cutoff is Ae- For pE (QCD), the IR cutoff is 
Ae, and the UV cutoff will be taken to be Jje (~ 27rT), i.e., the (MS) renormalization scale for gs(Ji^) appearing in pE- 



The same procedure gives for pM the foUowing expression: 

/^\ f^ 1 /^2 

1 



PM(r) = TmU ^ W 1 + :i3' |!L 



4 V2mEy 



+ 



2^2 

33 / 5e 



(47r)2 



Toe 



89 7r2 11 ^ 
24 6 6 



1 / 2 \ 3 

1 34 f 5e 



Toe 



(47r 

3 (-5) aW 
47r 2 ?7iE 



(q:m + "g) In ( j - 8 ac In 



2tot: 



+ /3] 



M 



(22) 



The expression in the brackets containing Q!m and ctQ is obtained from the term 8 au ln(/i/(2TOE)) in the expansion 
for pm/{TJ1~'^'^) of Ref. 0|, by replacing /I in part of the logarithm by the corresponding UV cutoff Ae of EQCD and 
in the rest of the logarithm by the IR cutoff Am in such a way as to guarantee that the sum pQ + pM is independent 
of Am- For this, relation H18|) for ^m must be inserted into expression (|19|l 



PG == TtoI 



(47r)3 



2 \ 3 

rriE 



Q 1 I ^M 

8 ae In -r^ 
6.91; 



(23) 



In expression (|22|l for pM, the values of constants au and /3m have been obtained in Refs. U and Q, respectively 

43 491 



aM = ' 

32 6144 



0.555017 



/3i 



'M 



-1.391512 . 



.(2) 



The last term in expansion H22|) is written in the convention Ag = were A^ is equal to [34 

1 

'24^2' 



X^^^^Tgti-p):r^{9-nf) + Oigt). 



(24) 
(25) 

(26) 



Ji) 



Here, Uf is the number of active quark flavors. In ref. |35l | a different convention is adopted yielding Ag — 
Tgfn//(4327r2) + 0(gf ) and AJ^^^ = Tgj{l - n//3)/(247r2) + ©(gf), but this then leads to the same result for pM- 

Adding expressions l|22|l and l|23|l. we obtain the sum pm+g as expansion in powers of the EQCD parameters gE 
and rriE 
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The other matching coefficients ttie and ^e can be expanded in powers of QCD coupling gg — gs{fJ.) (cf. 4], with 
e^^O) 

2 " 



m 



E 



T\l + nf/6)gUl+ (|^) P„(«/) + 2/3oln f ^j 



0(5.^ 



+ 0(.9^) 



.I^T,.^|l+(|^) [p>,) + 2,oln(^) 

where /?o = (11/4) (1 — 2n//33) is the first beta coefficient, and 

(0.612377 - 0.488058 nf - 0.0427979 nj) 



Pm{nf) 



(l + n//6) 
-0.387623- 0.423454 n/) 



(28) 
(29) 

(30) 
(31) 



The Ag /mE-term in expansion H27|l can be written with the help of the leading order relation (|26|) and relations 
(|28|l - (|29|l in several ways, for example in the following two ways: 

3 {-5)X'i^ _ 5 {9-nf) 
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rriE 
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The more conservative approach in the resummations based on EQCD expansion (|27|) should be to resum separately 

the expansion of terms which are powers of g'^/m^ and the expansion of terms which are powers of Ag /niE, because 
the two expansions represent probably two topologically different families of diagrams. The problem with the latter 
expansion is that we know only the leading term there, which is written in terms of the QCD coupling parameter 
gs(j£) in Eq. H32|l . and in terms of the first EQCD coupling parameter g^/rriE in Eq. I|33(l . The latter equation, in 
comparison to the former, represents some kind of (EQCD-)resummation of the Ag /rriE-family of terms. Since the 
parameters g^/niE and Ag /toe are both of the effective EQCD theory, the expression (|33|l may be considered as 

better, but less conservative. While Ag /ttt-e should be independent of the renormalization scale Jl = ^m (~ "t-e), we 
will later see that the leading order EQCD expression 133|) is significantly less /iM-dependent than the leading order 
QCD expression (|32ll . especially when the Pade resummations P[l/l](a(/iM)) are applied to expansions (|28|l - 129|) for 



and 



9e- 



The coupling parameter ^e and the Debye screening mass ttie are physical quantities, and thus they arc independent 
of the renormalization scale (/l). This independence is reflected in the coefficients of expansions H28|) - (|29|l in the terms 
proportional to PolnJI. Inserting expansions H28|l and (|29|l into series (|27|l . and using relation H26|l . leads to an 
expansion for the sum pm+g in terms of the QED coupling gs = gs (jl) 
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(34) 



where the constant K^ is 



K^ = 



12(l+n//6) 



(-99.9089- 35.1402 Tiy 



7.08145 nj) 



(35) 



The dependence on the factorization scale Am in the expansion H34|l disappeared as it should. Further, the quantity 
Pm+g is independent of the renormalization scale /I, and this is reficcted by the term proportional to /^oln/I in the 
coefficient at g^ jn Eq. (|34|l . 

We note that the only unknown dimcnsionless coefficient in (|34|l is Sq which, as argued before, is expected to be 

I'5g|^1. 

The knowledge of the expansion of the short-distance pressure p^ is less complete ~ this is an expansion in powers 
of (7^, and it is known only up to ~ g^ (cf. Ref. |3]). However, important parts of the coefficient at g^ in pE can be 
deduced from the requirement of /I-independence ofpE and of the AE-independence of pe+m+g (Ae is the factorization 
scale for the sum pE +Pm+g)- After some (tedious) algebra, we end up with the following expansion of pE up to ^ g^: 



Pe{T) =Pidcal(r) 



15 (1 + 5n;/12) „ 
4 (l + 21n//32)^ 
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(1 + nf/6)i3P„,{nf) + X3/6) + 18(aM + ^g) 



0{gt)\ ,(38) 



and the parameter k was introduced such that the NLO coefficient is made up of only two logarithmic terms propor- 
tional to Iii{Jl/2ttT) and \n{AE/KT) as shown above 



"^ K ) 135(l+n//6) 



(244.898 + 17.2419 n/ 



0.415029 n^) 



(39) 



For example, k ~ 1.0241, 1.47922 for Uf = 0,3, respectively. At order g^ in expansion (|38|l . only the dimensionless 
number S-^, is unknown - it is independent of the energy scales and of their ratios, just like i5g- We organized the 
coefficient at g^ in Eq. H38|) in the following way: the /I-dependent terms are written in powers of \n{'p/2TrT); the IR 
cutoff (AE)-dependent terms are written in terms of ln(AE/KT), because this combination absorbs all the ln(jl/2TrT)- 
independent terms in the coefficient at g^. It is reasonable to expect that the ln(7l/27rT')-independent terms at ~ g^ 
would also be absorbed to a large degree by a quantity proportional to the combination ln(AE/KT). That's why we 
expect the number Se to be small. A conservative expectation would be that \Se\ is not larger than the ln(AE/Kr)-term 
there: 



„(o) 



|fcr(AE)| < Se < +\k^"'{A 



(0), 



where 



fcf (Ae)^ 



36 



(1 + 5n,/12) '^ (§) (-T2^^(^ + "^/^) + ^^("- + "-) - l^' + n,/6)PUnf) 



(40) 



(41) 



While Se is independent of any scale, k^ will have only slight dependence on temperature T when we will take 
Ae = ^2TTTmE{T) (- ^y^T). This choice of Ae was taken also in Ref. ^. 



III. NUMERICAL ANALYSIS OF UNPHYSICAL DEPENDENCE ON SCALES 



Expansions (|34|l or H27() for pm+g, and H36|l - H38() for pE, in conjunction with expansions H28|l and (|29|l for m^ and 
(7e, will form the basis of our numerical analysis. The only unknown parameters are 6g and i5e which will be allowed 
to vary in the intervals (|21|) and H40|l . The numerical analysis will be performed in analogy with that of our previous 
work 21]. There, expansion corresponding to present Eq. (|34|l for pm+g was applied but contained only terms up to 
~ 5s ; expansion for i?|f" of present Eq. lP5|l contained only terms up to ^ gf; and expansion l(^ only the leading 
term. The resummations in Ref. J2lj | were performed with Pade approximants, or with simple evaluation of the 
truncated perturbation series (TPS). In the present work, the resummations will be performed with Fade, Borel-Pade 
(cf. Appendix EI), or TPS. 

As argued in the Introduction and in Ref. |2l| . the MS renormalization scale Jl should be chosen accordingly in 
different regimes for the resummation of the different quantities (jSU, (|38|l . and H28(l - H29(l . The scale regime for JI 
should be of the order of a typical physical scale that corresponds to the quantity to be resummed. Therefore, the 
renormalization scale /I = /^e hi the short-distance quantity 138|) is /ie ^ 2t:T. For the long-distance EQCD quantities 
(P5|l - (j^ . the relevant scale is /I = /"m such that ^m ~ "m-E (~ gsT). For the long-distance quantity pm+g, Eqs. l|H^ 
and H27|l . the relevant scale JI = /im+g should be somewhere between '^ gsT and g^T; we will take it ^ gsT, i.e., 
JI = /ZM '^ gsT. Unless otherwise stated, we will take /iM = IJ^m = iTiE^l^-m)- For the factorization scale Ae we take, 
unless otherwise stated, the geometric mean between the hard scale 2ttT and the EQCD scale toe: Ae = [27rTTO,E]^/^ 

(~ 5s T). For to| and 5I we will take, unless otherwise stated, the P[l/l](a) Pade approximant (PA) with respect 



10 



to a (cf. Appendix IX|) . where a = a{ji) = {g s{ji) / IttY' and the scale /I (= /im) is adjusted so that /I = toe- We will 
see below that P[l/l](a) is a very reasonable resummation for m\. 

Further, if not stated otherwise, we will take for the number of active (and massless) quark flavors Uf — 3. For the 
QCD coupling parameter asijl) = g^ (p) / {An) we take the reference value as(TO^, MS) = 0.334 which is approximately 
the value extracted from the hadronic t decay data |3a. l37| . We work in the MS scheme and use for the /3 function 
P[2/3](a) Fade approximant (PA) (a = Ug/n), unless otherwise stated. This approximant, as shown in Refs. [l7l |. 
represents a reasonable (quasi) analytic continuation of the TPS of /%jg(a) into the regime of large a [with JI down to 
as{iJ.) « 1.0] where the TPS is practically inapphcable. 
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FIG. 1; (a) The Debye screening mass ms, and (b) the EQCD couphng parameter g|,, as functions of the renormalization scale /im, 
when T= 1 GeV. The upper of the LO TPS (NLO TPS) twin curves has a(/i^) evolved by the one-loop (two-loop) RGE from a{m'^). 
The other curves have a{fi^) evolved by the four-loop PA [2/3] beta function. 



First, the resummation of the expansion of the squared Debye screening mass to| (|28(I and of the EQCD 
coupling parameter g| l|29|l . respectively, are performed. Both expansions are next-to-leading order (NLO) in 
a{fim) = [5s(Mm)/27r]^. Therefore, the diagonal P[l/l](a) can be constructed and should be a good candidate. 
Figs, n^ a) and PT^b) show the results, as a function of the renormalization scale fim (~ itie)- The values of toe in 
Fig. n^a) are obtained by the corresponding resummation of expansion H28|l for m^ and then taking the square root. 
In addition to P[l/l]{a{fim)), also the effective charge method result (ECH) |3a. l39l l4(]|| . as well as the TPS results, 
are presented. The ECH result is the NLO TPS at a specific value of the renormalization scale, so it is automatically 
independent of fim- We see from Figs, ^^a) andQIb) that PA's P[l/l](a) for r7i| and g| are very good approxima- 
tions as they almost eliminate /J,m-dependence contrary to the TPS-expressions which show (both for LO and NLO) 
a very strong unphysical /^-dependence. This is in accordance with Ref. [l3| where it is argued that any diagonal 
P[n/r7,](a(/i)) of a QCD observable has significantly reduced /i-dependence, i.e., it is /^-independent in the large-/9o 
limit.* As mentioned before, we will choose the mass ttie by the condition (to|)^1^/^1 



square root of this value as m^ ' (T) 



/(^, and we will denote the 



=(r)) 



p[i/i] 



= A*„ 



,(0) 



(T) 



(42) 



At T = 1 GcV (and Uf = 3), this value is to^°^ « 1.923 GeV.'^ 

In Figs. 121 we present the short-distance pressure pe as a function of the respective renormalization scale /xe - 
various resummed expressions based on the perturbation series (|5^ - (|55|) in powers of a{^E) = [5s(ME)/27r]^, for T = 1 
GeV and (5e = 0. The naming of each resummation refers to the use of the corresponding approximant for i?|f" in 
Pb of Eq. (|36|l as a function of a(/iE) = [5s(ME)/27r]^. For example, P[2/l] means that we apply to the perturbation 
series (|38|) in powers of a(/iE) the Fade approximant P[2/l](a(/iE)); BP[l/2] means that we apply P[l/2] to the 



** There exists an extension of the diagonal P A's such that it is completely /^-independent llfil , or fi- and scheme-independent llTl . 
^ Fig.0^a) was presented in our previous work |21|. but there the curves are slightly lower (by about 0.02 GeV) due to inadvertent omission 
of the factor (1 + nf/6)~^ appearing in our Eq. 1301 for Pm- 
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FIG. 2: The short-distance pressure pg, at T = 1 GeV, as a function of the corresponding renormahzation scale fi^ - various approxiniants 
are applied to the perturbation expansion 1381 in powers of a(fiE) = [9s(/JE)/27r] , and <5e = was taken. Figure (b) has a finer vertical 
scale and includes some approxiniants not shown in Fig. (a). 



Borel transform of expansfon H38|l - cf. Appendix 1X1 The range of /ie is around 2ttT, i.e., the order of the relevant 
physical modes contributing to pe- Further, the IR cutoff Ae is taken fixed: Ae = [ZnTm^' (T)]-^/^ . The results are 
normalized by pidcai = (87r^/45)(l + 21n//32), i.e., the expression for the pressure when T ^ oo [cf. Eq. (|36|l ]. From 
Fig. I3a) we see that the TPS's are not acceptable, due to too strong /iE-dependence. For the same reason, as seen 
from Fig. |2Ib), the approximants P[l/2], P[2/l], and Borel-Pade BP[2/1] are not acceptable. The only acceptable 
candidates are P[l/1], and BP[l/2]. We note that P[l/1] is based explicitly only on the NLO series (up to ~ a^) of 
Eq. H38|l . Nonetheless, we will include later P[l/1] of i?|f" in some of our results, because it is close to BP[l/2] and 
is very stable under variation of /xe (NL ECH also gives values close to P[l/1] and BP[l/2]). Although the results 
are presented only for the case T = 1 GeV and ^e = 0, the behaviour and the conclusions remain the same for other 
temperatures, and for other values of Sb in the interval (|40|l . 



o 
+ 

Q. 





, pn /1 1 _ _ _ 1 


0.4 








- 






BP[1/1I ■ 




0.2 


(a) 




- 






P[0/31 











- 




~TPS 




n, = 3, 5q=0 
T = 1 GeV 


0.2 


P 
p 
p 


0/3 
2/1 
1/1 


77^..^ 




Ae = 3.476 GeV- 


0.4 


■ BP 
BP 


2/1 
1/1 






p[2/i] :" 

,---'-'"-' -,--BP[2/1] 



1.5 



2.5 
V^M [GeV] 



3.5 



^ 

^ 



0.2 

0.15 h 

0.1 

0.05 



-0.05 - 

-0.1 - 

-0.15 

-0.2 
1 




n, = 3, 8q=0 



TT1/2] 



T = 1 GeV 



= 3.476 GeV 



P 

P 

BP 

BP 

BP 



1.5 



2.5 
Hm [GeV] 



3.5 



FIG. 3: The long-distance pressure pM-i-Gi at T = 1 GeV, as a function of the corresponding renormalization scale ^m — various 
approximants are applied to the perturbation expansion Il34i in powers of gsifJ-m)^ a^nd 5q = was taken. Figure (b) has a finer vertical 
scale and includes some approximants not shown in Fig. (a). 



Similar analysis is now performed for the long-distance pressure pm+g- Pade and Borel-Pade resummations are first 
applied to expansion 1)34(1 forpM-i-G (oc Pm+g/it^e) which starts with one and is in powers of (7s(mm) ~ cf. Appendix IXIfor 
more details. In Figs.|31we present the results of various resummations as a function of the respective renormalization 
scale /iM, for T — 1 GeV and Sg = 0. The factor to^ in the first line of Eq. H34|l was taken with P[l/1] for m'^dim), with 

the renormalization scale fim = Mm- The UV cutoff Ae was taken fixed according to the formula: Ae = {2iTTm^ )^/^. 
The TPS and some of the other resummations show significant ^M-dependence. Further, BP[0/2] and P[0/2] are of 
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FIG. 4: Same as in Figs.|21 but now Pm-i-G is based on expansion I^YI in g|,/mE- Details are given in the text. 



terms of expansion H34|l 



^ in Pm+g)- We decide at this stage to consider only the N^LO resummations: Pjl/2], P[0/3], BP[l/2], 
The approxiniant P[0/2](gs) for Pm(+g) "^^^ considered in our previous work |21| . and we include it 



lower order (NNLO) and do not take explicitly into account the numerically important ~ g 

of Pm+g (~ . 

and BP[0/3] 

later in Figs.[7|and|51in the presentation oi p{T). The conclusions of this paragraph survive when other values of T 

and ^G are used. 

We alternatively apply our resummation procedure to the EQCD expansion H27(l instead of H34|l for pm+g- The 
results are presented in Figs. 01 The values of g^ and to| in (|27|) are chosen as P[l/l](a(/i„i)), due to their /.tm-stability 
(Mm, = Mm taken) as seen in Figs.^ and then Pade or Borel-Pade are applied to expansion 127|l in powers of the EQCD 
parameter g^/rriE {^^ gs), without the Ag /?7iE-term which is then added separately as the leading order QCD term 
ex 5s (mm), Eq. l|^ . When Pade or Borel-Pade based on TPS terms of order lower than ^ g^ for pm+g of Eq. l(77|) 
are applied - such as P[l/1], BP[1/1], P[0/2] or BP[0/2] to TPS ^ oi pm+g " the A^^V™E-term is not included 
and not added as it represents a term ~ g^ of pM-fG (~ 9s to pm+g)- The TPS curve in Fig.^a) was obtained by 
evaluating expansion (|27|l of pm+g directly as EQCD TPS, with the values for g^ and to| taken as simple (NLO) 
TPS's (|28|I - H29|I at Jl — /im- The conclusions from Figs. 01 are the same as in the procedure leading to Figs. 01 the 
N'^LO resummations P[l/2], P[0/3], BP[l/2], and BP[0/3] all remain acceptable candidates at this stage, and this 
conclusion turns out to be independent of T and Sq . 
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FIG. 5: The total pressure PQCD = PE -I- PM+G, at T = 1 GeV, as a function of the corresponding factorization scale Ae — various 
approximants are applied separately to the perturbation expansions i?g^"(a(^E)) I38i and Pm-i-gCSsImm)) 1341 . Details are given in the 
text. Figure (b) has a finer vertical scale. 



Another necessary condition for an acceptable resummation is that the spurious AE-dependence in the total sum pe+ 
Pm+g be weak. In Figs. [SI we present the results on AE-dependence, at T = 1 GeV, combining various aforementioned 
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FIG. 6: Same as in Figs.|3 but now Pm+G is based on expansion I2YI in g^/m^. Further details are given in the text. 



resummations for pE and pm+g that were found acceptable so far. Here, toe — rri^ was fixed in the aforementioned 

way H42|) . he — ZttT, and /im — fJ-m = Wg . Now, Ae 7^ ^/MeMm, but varies between /um (~ 1-9 GeV) and /iE 
(« 6.3 GeV). The two unknown parameters were chosen to be 6e = 5q — 0. We see that the curves with P[0/3] and 
B[0/3] used for pm+g have acceptably weak As-dependence, while those with P[l/2] and BP[l/2] used for pm+g have 
unacceptably strong AE-dependence. Thus, we are left with just two types of resummations which fulfill the necessary 
conditions of weak /iE, /iM and Ae dependence: 

1. P[0/3](g,(AiM)) for pm+g; and either ¥[l / l]{a{iiE)) or BP[l/2](a(/iE)) for i?|f"(a(ME)); 

2. BP[0/3](g,(/iM)) for pm+g; and either P[l/l](a(/iE)) or BP[l/2](a(A*E)) for i?g'"(a(A*E)).^ 

These conclusions do not change when the values of T, 5e and Sq are changed. We note that P[l/l](a(/iE)) is of lower 
order and therefore does not use the term ~ a^ in expansion (|38|l . 

We can repeat the same analysis, but using expansion l|27|l for pm+g in powers of (.^e/to-e) instead of expansion 

()34|l in powers of gsdiu)- The values of g^ and jn^ are taken as P[l / l]{a{fiM)) , with /im ~ Hm — 
are given in Figs. El they are similar to those of Figs.|Sl and the conclusions are the same. 



,(0) 



The results 



IV. NUJVIERICAL RESULTS AS A FUNCTION OF TEMPERATURE 

We will first concentrate on the first family of resummations, i.e., those with P[0/3] forpM+G- We present the results 
for these resummations as a function of temperature T in Fig. Ufa), for Se = 5g = 0, and pm+g is based on expansion 
()34|l in powers of gsipu)- In addition to these resummations, we include also the same type of resummations where 
the pq part is excluded (pe +Pm), where we use in pM for the IR cutoff: Am — (to.^ )^/Ae. We use iie — 27rT; 



Mm = fJ-n 



Toe 



to 



(0) 



(T); Ae — y^/IeA^m- We can see in Fig. (Tfa) that the presence of pc, at least for 6q ~ 0, 
decreases the value of the total pressure somewhat. For comparison, we also include the result of resummation of 
P[l/1] for i?|™ and P[0/2] for pu (cf. Ref. |2j|), i.e., the case where the terms ~ g^ in pm+g (terms ~ gf in pm+g) 



are not explicitly accounted for (and neither are the terms 



5. 



in Pe). Fig. El^a) shows one interesting feature: 



when the terms ~ g^ in pm+g are explicitly accounted for in the resummation, the sign of the curvature becomes 
negative in the entire temperature interval - i.e., at least there where the resummation is applicable: T > 0.3 GeV. 
Thus the curvature at low temperatures has now the same sign as suggested by the known relation p/pidcai << 1 
at r ~ Tc « 0.2 GeV (see later in this Section). For T < 0.3 GeV our resummations cannot be applied any more, 
because in that case iJ,m{= Mm) < 0.81 GeV, but the P[2/3](a) beta function does not allow running below such values 
[a{iim) blows up]. When we vary the values of Sq and 6e in the intervals 1^1)1 and iPHjI . respectively, the predictions 
do not change much. This is presented in Fig.[7|[b). 



NL ECH could be used instead of P[l/1] or BP[l/2] for i?5=f", but the results are similar in all three cases, at any temperature. 
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FIG. 7: (a) The total pressure p (normalized by Pidcal) ^s a function of temperature T, with 5^ = Sq = - various resummations 
are applied: for R^'^{a{fi^)) the approximants P[l/l](a(/.tE)) and BP[l/2](a(/iE)); for pm+g(9s(mm)) the approximant P[0/3](gs(/.tM))- 
Shown are also the analogous results when pQ is excluded. In addition, a resummation which does not account for the ~ gf terms in 
Pm+G is included (dotted line), (b) Variation of a specific Pade/Borel-Pade resummation when the unknown parameters Sq and 5e are 
varied. The full curves in (a) and (b) are the same curves. Further explanations are given in the text. 
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FIG. 8: Same as in Fig. 171 but now pM+G is based on expansion I27i in g|,/mE - using the resummed g^ and m|, as P[l/l](a(^M))- 



When we do not base the resummations of pm+g on expansion (|34|l in powers of gsifim), but rather on expansion 
(|?7jl in powers of the EQCD parameter g^/niE [using the resummed g^ and m| as P[l/l](a(/^M)), and adding the 
•^E /"^E-term separately as the leading order QCD term ex g'^ifJ-u), Eq- ISSl], remarkable changes occur for the results 
of p/pidcai as a function of temperature. The obtained results, analogous to those of Figs.Cfa) andEl^b), are presented 
in Figs.lHfa) and|Htb), respectively. The "low order" dotted curve in Fig. |Hla) (the uppermost), which does not use 
information on ^ g^ terms (and thus does not include the Ag -term), still does not change the curvature in the low 
temperature regime. On the other hand, "higher order" Fade resummations of pm+Gi which use information on ^ g^ 
terms in pm+g {^ 9^ in Pm+g) and include, added separately, the Ag /rriE-term in pm+g as the leading order QCD 
term ex 5^ (/iivi) [Eq. (j32(l ]. result in a pronounced negative curvature and a rapid fall of p/pidoai when the temperature 
falls down toward the critical values Tc ~ 0.2 GeV. This behavior is qualitatively correct because we know that it must 
be p/pidcai << 1 at r « Tc (see later). In that respect, this phenomenon indicates that the resummations of pm+g 
based on expansion H27|l in powers of the effective EQCD theory parameter g'^/m-^ are more reliable than those based 
on expansion H34|l in powers of gs- The low-T results of Figs. |H1 turn out to be closer to those of lattice calculations, 
the feature which will be discussed and presented in more detail in Section IVII Another very positive feature can be 
read off from Fig. Etb) : variation of the Fade and Borel-Pade resummed predictions, when the unknown parameters 
5g and <5e are varied in the generously wide intervals H21|l and H4()(l , is weak. 

Now we turn to the second family of resummations, i.e., those with BP[0/3] for pm+g instead of P [0/3]. We present 
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these results, as a function of temperature T, with values of i5e and ^g varied in the intervals H40() and HU, in 
Figs-inia) andlHI^b). In Fig. |5fa), BP[0/3] was appHed to expansion (|Sl|l for pm+g; in Fig- Elb) to expansion |[77|) 
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FIG. 9: (a) The total pressure p (normalized by Pidcal) ^s a function of temperature T, for various values of <5e and <5q, when BP[0/3] 
is applied to expansion I34i for pM+G (and BP[l/2] to i?^^"); (b) Same as in (a), but BP[0/3] for Pm+G is based on expansion I27i in 



g^/rriE, using the resummed g|, and m|, as P[l/l](a(/iM))- 



for pm+g without the Ag /rreE-term, the latter was taken as the leading order ex 5s (/^m) of QCD expansion (|32|l and 
added separately (only in the "higher order" P[0/3] and BP[0/3] cases). Further, the other parameters and procedures 
are the same as in Fisfs. 17181 For comparison, the lower order counterpart, i.e., with BP[0/2] for pm+g (and P[l/1] 
for i?g*") is included in these Figures. In addition, for comparison, the corresponding curves when P[0/3] is applied 
for pm+g, and at lower order P[0/2], are included whenever visible^ - see also Figs. 17181 We observe from Figs. Inland 
I7l8l that the choice BP[0/3] for pm+g gives p/pidcai > 1 for most of the low temperatures T, and the choice P[0/3] for 
Pm+g gives p/pidoai < 1- 

We recall that BP[0/3] for pm+g was equally acceptable as P[0/3], when only the demand for weak renormalization 
and factorization scale dependence is used as a criterion, Figs. I3I6I However, there are several indications (not using 
the approximate knowledge of the low-temperature curves of p/pideai from lattice calculations) that the choice BP[0/3] 
is less acceptable than P[0/3]. Some of them are due to general physical considerations, others are connected with 
the specific numerical approximation techniques applied here. 

Physical considerations provide at least two arguments for expecting p/pideai < 1 and, consequently, for favoring 
the choice P[0/3]: 

1. From the physical point of view one would expect that the pressure of a (relativistic) quark-gluon gas gets 
lowered relative to the free particle case once the interaction is switched on, simply because we expect that the 
behavior of interacting massless quarks and gluons is approximately described by (almost free) massive quasi- 
particles, the mass stemming from Debye screening. Such a behavior is manifest not only within a non-relativistic 
electromagnetic plasma (calculated according to Debye-Hiickel) , but shows up also in specific model calculations 
for a relativistic plasma ^ ,3 ■ Therefore, we expect that p/pideai < 1 for T close to critical temperatures Tc 
(w 0.15 -0.25 GeV). 

2. The same qualitative behavior is inferred from a thermodynamic consideration [41|: We know that the pressure 
(considered as a function of T) remains continuous at the phase transiton point T = T^. But below T^ the 
system is a hadron (mostly pion) gas, which to a good approximation can be described by an ideal pion gas. 
The corresponding pressure is 

TT^ 3 Stt^ 



It is natural to regard P[0/2] and P[0/3] (for Pm+g) ^-^ part of a sequence of approximations, and BP[0/2] and BP[0/3] as part of 
another sequence of approximations — e.g., Ref. Illl and Figs.l^Ja) andlWbV 
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TABLE I: The exact and predicted coefficients Tj of the expansion 1341 for pm+g in powers of gs{^J.M)- The scales are fixed as 
usual: MM = Mm = "^e' of Eq. ^, Ae = (27rrm^°V''^. The temperatures are T = 1 GeV and 0.5 GeV; Uf = 3. The resuhs 



''3(<5g) are given for the values 6g 


= ± 5 Eq. jnj 








T [GeV] n (exact) 


r2 (exact) 


r:i{5G) (exact) 


rs (pred. P[0/2]) 


r-3 (pred. BP[0/2]) 


1.0 -0.3794 
0.5 -0.3374 


-0.4654 
-0.4654 


-0.0643 ±0.1111 
-0.0852 ±0.1111 


0.4077 
0.3524 


1.387 
1.172 



TABLE II: The exact and predicted coefficients r^ of the expansion 1)27^ for pm+g in powers of the EQCD parameter g^/niE 



without the A^ ' 


-term. All the other parameters as in Table Q 






T [GeV] 


rf (exact) r^«(exact) rf{SG) (exact) 


rf (pred. P[0/2]) 


rf (pred. BP[0/2]) 


1.0 
0.5 


-0.4647 -0.6980 -0.1645 ± 0.2041 
-0.4132 -0.6980 -0.1902 ± 0.2041 


0.7490 
0.6474 


2.548 
2.154 



which is much smaller than the pressure for an ideal gas of quarks (with Uf flavours) and gluons, Eq. H37|l . 
Consequently, the true pressure of the interacting plasma at temperature T > Tc (but close to Tc) must be 
much smaller than the ideal gas value at the same temperature, and we again obtain p/pidcsd < 1 for T ^Tc. 

From Fig. |5Ja) it is seen that p/pidcai > 1 above Tc and even seems to increase when T —> Tc if the approximant 
BP[0/3] is applied to expansion jM|) of pm+g (5e ^°t being used/resummed), whereas P[0/3] is in accordance with 
the aforementioned expectation that p/pidced < 1 [Fig-I3a)]- On the other hand, when expansion (|27|l is used as the 
basis for rcsummation (with g^ resummed) , p/pidced falls below one for temperatures very close to Tc even in the case 
BP [0/3] (but g^ resummed) , but it is still larger than 1 for a considerable temperature region above but not far from 
Tc [Fig.infb)]. This indicates that using expansion (|27|) instead of H34|) gives in general more realistic results. 
Numerical considerations provide at least two other arguments in favor of the choice P[0/3] for pm+g'- 

1. If BP[0/3] for Pm+g gave results closer to the true values of pm-i-g than P[0/3], then Fig. I^b) would suggest 
that the lower order counterpart P[0/2] (to P[0/3]) for pm+g gives results which lie closer to the true values of 
Pm+g than those from P[0/3] for T > 0.5 GeV - a situation that has to be regarded as unlikely.* 

2. The 0{g^) terms of pm+g predicted by re-expansion of BP[0/2] are clearly worse than those predicted by 
P[0/2], cf. Tables 1 and 2. This suggests for the respective higher order approximants BP[0/3] and P[0/3] the 
corresponding hierarchy of reliability. 

For all these reasons, we will regard as the acceptable resummation at the 0{gg) level to be the one using P[0/3] for 
Pm+g, and BP[l/2] for i?|f". Further, as mentioned before, the resummation oipM+G appears to be more consistent 
with the physical expectation of p/pidesd < 1 s^t T —^ Tc when it is based on expansion H27I) of pm+g in powers of the 
EQCD parameter g^/ruE (with the Ag /wE-term added separately). 

Up until now, all the results presented were for nj = 3. In Figs. ^| we present comparison of (Borel-)Pade 
predictions in the cases Uf — and Uf — 3: in Fig. llOf a) when pm+g is based on expansion (|34|l in powers of gs; and 
in Fig. llOf b) when pm+g is based on expansion (|27|l in powers of g^/niE, with the aforementioned treatment of the 

Ag /TOE-term. In the latter Figure we also included two curves for the case when the MS /3 function is taken as TPS 
instead of P[2/3](a(/i)). We see that the negative curvature at low T survives also in the nj = case. Further, the 
results for n/ = 0, 3 depend only very little on the type of the resummation used for the /3-function, even at low T. 

We present the main results for p/pideai as a function of temperature, which are given also in Figs.|HIb) and lTUT b). 
in a detailed form in the low-temperature regime in Figs. llll In the latter Figures, we present further the variation of 
the curves when the renormalization scales /xe, and hm = Mm are varied around their central values 27rT and m^ (T) 
of Eq. (|^ by factors 1.5 and 1/1.5.^ The variation of /Ue changes the curves insignificantly. 



* Note that the lower order "counterpart" P[l/1] (to BP[l/2]) for R|?" gives, at any temperature, values of pE very similar to those of 

BP[l/2] - cf. Figs.|3a) andlUa). 
^ The factorization scale is Ae = {2nTmE)^^^ where toe is the square root of P[l/1] of r?i|,(Mm). Ae changes only little with the variation 

of the scale fiyi = ^m- 
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FIG. 10: (a) Analogous to Fig. |7Ia), but with nj = for comparison; (b) as in Fig. (a), but now Pm+G is based on expansion 1271 in 
g^/rriE — using the resummed g|, and m,|, as P[l/l](a(^M))i i.e., analogous to Fig.lHJa) but with nj = for comparison. 
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FIG. 11: (a) The pressure p, normalized by Pidcal, a-s a function of temperature T, for the most acceptable resummation approach: P[0/3] 
for pm+G of expansion I27i . and BP[l/2] for flji?'", for (5e = 5g = 0, and nj = 3,0. The variation of the renormalization scales fiE and 
Mm = Mm is by factor 1.5 and 1/1.5 around the values 2nT and m^'{T) of Eq. 142 i . respectively, (b) Same as in Fig. (a), but now the 
unknown paremeters 5e and 6q are varied, according to Eqs. 1401 and I21i . respectively. 



However, the variation of the lower scale /im = Mm influences significantly the results for pm+g and thus p/pide&h 
at T < 1 GeV, as seen in Figs. ^J This strong variation is a reflection of at least two aspects present at T < 1 

GeV: (a) the scale /im = /^m ^ "t-e (^) falls down to < 1 GeV, a region where the perturbative approach and RGE 

running eventually break down; (b) the hierarchy of scales (SttT) 3> m^ (T) gets very narrowed down. The first 

aspect is reflected in the strong instabihty of the leading order QCD term oc g^ipm) [Eq- (|^ ] for the A^ /mE-term 
in expansion l|27|l under the variation of /im (= (Xm) at low T's. This term is not included in the resummation, as 
emphasized earlier. The variation of this term is a major source of the appreciable variation of the p/pideai-curve 
at low T's in Fig. Ilir a). The variation of jjm = fJ-m downwards to 771^/1.5 = 0.73 GeV is not allowed because 
the coupling parameter gs{fJ-m) blows up at such low scales by the renormalization group equation when the beta 
function is continued into the strong coupling region by the Fade [2/3] (as is the case in our numerical results). 

The renormalization scales /iM = M?ti = m^ /1.5 correspond to the lower curves in Fig. 1111 and they were drawn 
down to such temperatures where the corresponding coupling constant a(fj,„i) blew up. For the other renormalization 
scales, the curves were drawn down to approximately such temperatures where the physically motivated condition 
(|42(l cannot be fulfilled any more. 

We can apply for the Ag /?7iE-term in expansion H27() of pm+g the leading term of expansion H33|l in powers of the 
first EQCD parameter g'^/m^ (with P[l/1] for g^ and P[l/1] for 77i|, at renormalization scale /iivi), instead of the 
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FIG. 12; Analogous to Fig. Illl a). but now the results are included when the Ajj /m,E-term in expansion I27i is written as proportional 
to (gg/niE)'' [Eq. I33i ]. and not gj [Eq. I32i 1. The variation under the aforementioned changes of the scale /^m = ^im are included; 
Se=Sg= 0. 

leading term of expansion (|32|l in powers of the QCD parameter 5s (mm)- The X^ /niE-ieim is again not included in 
the [0/3]-Pade resummation, but is added separately. These results are shown in Fig. EJ along with the results of 
Fig. 111! Only the variation of the curves under the aforementioned variation of the /zm = Mm scales is presented. The 
variation of the two types of curves, at any given temperature T, under the changes of /xe, Se and Sq, is identical, thus 
very weak (shown in Figs.lll|l. We see from Fig. ^1 that the variation of the curves under the changes of fiM = /im is 
now significantly weaker. This has primarily to do with the significantly weaker /iM-dependence of the EQCD term 
(5|/toe)'^ of Eq. (j2ni), in comparison to the /iM-dependence of the the QCD term g^ifJ-u) of Eq. (|S^ . 



V. EVALUATIONS OF TRUNCATED PERTURBATION SERIES 



Now we investigate how the results change if, instead of Fade or Borel-Pade, simple TPS evaluations are applied to 
all expansions: to expansion l|^ for i?g"\ to l(^ and l(^ for to| and g^, and to l|5^ or (|77|l for pm+g- The results 
for the pressure, as a function of temperature, for 



(|5^ and (|?7jl are used for pm+g, respectively. 



= 3, are presented in Figs- ll^f a) and ll3f b). when expansions 

2\(TPS)/'.. ^^l/2 



Uf 

We enforce here: 



fJ-M = fJ-m'i and either /i„ 



[ir 



Ht^m)] 



or M„j == [(™|)^'^/^'(Mm)] [the latter is m^^' (T) of Eq. ^]. The curves with (NL)TPS m| are continued down 
to such temperatures where the aforementioned condition fim — fnE cannot be enforced any more, indicating that 
this (perturbative) TPS method becomes inapplicable below such temperatures T « 2 GeV. The curves with P[l/1] 
77i| in Figs. 1131 can in principle be shown, as earlier, for temperatures down to T « 0.3 GeV; however, their values 

(for |(5e| < 1^2 I) fall drastically: p/pidca.\ < 0.45 already at T w 1 GeV (which, incidentally, is far lower than the 
lattice results), and p < already at T « 0.5 GeV. These curves, with P[l/l]-resummed to| and g| are presented 
in Figs. [T^ in more detail for the low temperatures, where now Uf = 3 and Uf — 0. The A^ /?TiE-term in expansion 
(|27|l is added as the term c>c ggifJ-u) [Eq- H32I) ]. In Fig. I14f a). the variation of these TPS curves under the changes 
of the two renormalization scales fiE and /im — IJ-m, as explained for Fig. Illf a). is presented. In Fig. I14f b). the 
variation under the changes of the parameters (5e and ^g is presented. In Fig. I14f a) we see that the TPS curves vary 
more strongly under the changes of the high-energy scale fiB than under those of the low-energy scale fiu — f^m- 
This result, at first sight paradoxical, occurs mainly because we used for pm+g the EQGD expansion H27|l . with the 
dominant part coming from powers of the EQCD parameter g^/m-E, where the /^m-dependence of this parameter is 
rather weak because we used P[l/1] (a(/im)) for g^ and for rn^. Comparing with the corresponding Pade-Borel-Pade 
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(P+BP) resummed curves of Figs. ^2 we see that the latter have much weaker dependence on the parameters Sb and 
Sq and on the high-energy renormahzation scale /^Ej while the dependence on the low energy renormalization scale 
Mm = Mm (~ fnE ^ QsT) is reduced by P+BP resummation only by a factor of 2-3. For example, at T = 0.7 GeV and 
Uf = 3, this variation is about 0.18 and 0.06 GeV in the TPS and P-I-BP cases, respectively. This has largely to do 

with the Ag /"^E-term which is not resummed in the P-fBP case, but is also taken as ex 5s(mm) [Eq. (|H2)) ] and added 
separately. On the other hand, taking that term as oc (gl/TTiE)^ [Eq. H33|l ]. and adding it after the resummation, 
significantly weakened the /iM-dependence at low temperatures, as was seen in Fig. 1121 

On all these grounds, the Pade and Borel-Pade results of Figs. II II and 1121 are likely to give more realistic results at 
low T ^ 1 GeV than the corresponding TPS results. 
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FIG. 13: (a) The total pressure as a function of temperature T when TPS evaluation is employed for iJS?" and Pm+Gi instead of the 
Pade and Borel-Pade of Fig.|3 (b) Same as in (a), but now pM+G is evaluated as expansion 1271 in g^/m-E,. Further explanations are 
given in the text. 
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FIG. 14: (a) The low-temperature total pressure (for Uj = 3, 0) when TPS evaluation is employed for i?.^^" and Pm+Gi but for 

g|, Pade P[l/1]. The variation of the renormalization scales /xe and fim = A^M is by factor 1.5 and 1/1.5 around the values 2ttT and 

m^'{T) of Eq. 142 i . respectively, (b) Same as in Fig. (a), but now the unknown paremeters 5^ and 5q are varied, according to Eqs. l4Ui 
and I21i . respectively. 



VI. COMPARISONS WITH OTHER APPROACHES, AND CONCLUSIONS 



An approach different from our resummation is to set all (MS) renormalization scales equal: ^ = fi-^ = Mm = Mm = 
Ae, and then evaluate ttie, pe and pm+g (=^ Pm+g)- This was the approach of Ref. 42], and the evaluation in 
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FIG. 15: (a) The total pressure (with ny = 3) as a function of the temperature T when TPS evaluation is employed for i?|?", m,|,, and 
Pm+G of 1341 . but with common scales used: /xg = A'm = Mm = ^E = l-797rT. The unknown parameters Sq and <5e are varied in the 
interval 1211 and 1401 . respectively, for nj = 3. (b) same as Fig. (a), but using for pM+G expansion 1271 in powers of g^/ms,. Further 
explanations are given in the text. 



^ 

^ 



0.9 
0.8 

0.7 
0.6 
0.5 
0.4 



0.3 





(a) _-^^^-" -' 




n,=0, |x=1.79jtT 




- TPS's: Rg", Pm+g ,/■'■ 
m|isNLTPS 




g|isNLTPS / 
■ S, = S(;=0 / 

/ y 


// / 


A ■' 
A ■' 
A ■ 
/ / 

//■■■ , 


/ 



0.3 



0.5 



^ 

^ 



1 
T [GeV] 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 



n,=3, 6e=5q=0 
nf-3, var.of 5q= ±-5-"" 
n,=3, var.of,'6'^= ± kf 
n,=0,,&E=6Q=0 
rif^t), var.of 5q= ±5 
■''n,=0, var.of 5p= ± k!?' 



/\^'E= \^m= V'U= f^e^ 1-79 ItT 
TPS's: Rg", Pm+g 



niE is NL TPS 
g| is NL TPS 



0.3 



0.5 




FIG. 16: (a) The low-temperature total pressure (for n^ = 3,0) when TPS evaluation is employed for fljjf", fn'^, 9^, and pM+G of 
Eq. 1271 . but with common scales used: /^e = Mm = Mm = Ae = l.TdnT. Variation of the scale is by factor 1.5: Mmax = 1.51.797rT, 
Mmin = (1/1-5)1. TQttT. (b) The unknown parameters i5g and (5e are varied in the interval 1211 and 1401 . respectively, for Uf = 3,0. 



Ref. 3 was similar as well. In both of these references, truncated perturbation series (TPS) evaluations were applied. 
In Figs. llSt a) and llSf b'). we present the results of such type of resummation as a function of temperature, for Uf = 3, 
when pm+g is evaluated as TPS (|34|l and TPS (|27|l . respectively. The common scale was chosen to be p = 1.79ttT and 
the expansions (|38|l . (|28|) . H29() . (|34|l [or (|27|l ] were evaluated as TPS. This should correspond roughly to the method 
leading to the the dash-dotted curve in Fig. 5 of Ref. [ij for Uf — 0. We can see from Figs. E| that this approach 
gives, at low T < 10 GeV, results similar to the dotted TPS curves of Figs. IT!^ where P[l/1] was used for m^ and 
different renormalization scales were used for TPS pe and TPS pm+g- However, in contrast to Ref. |42|, most of our 
TPS curves in Figs. 1131 and 1151 fall down faster then theirs when temperature decreases to T < 1 GeV [T/T^ ^ 5), 
and fall down faster than the lattice curves. Only the upper Uf — Z curve in Fig.^Jb), corresponding to the choice 

5'E — ~\kA\ (and 5q = 0), is marginally compatible with the lattice results. 

In Figs. llfif a'l and (b) we present, at low temperatures, and for Uf — 3 and 0, the variation of these curves when the 
common renormalizatio scale ^ changes around 1.797rT (from I.TOttT/I.S to 1.5 x 1.797rT), and when the parameters 
(5e and 5q change, respectively (in analogy with Figs. 1141 and lll|) . We see that the variation of the curves is then 
similarly strong as for the TPS curves of Figs. E] 

On the other hand, Pade and Borel-Pade resummations. Figs. 1811 II (falling down as well at low decreasing T) are 
much less dependent on the unknown parameters 5q and (5e, as can be seen by comparing with the TPS curves 
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FIG. 17: Summary of our results for nj = 3 and Uf = 0, when the unknown parameters <5e and 5q are varied. The resummation of 
Pm+G was performed on the basis of expansion I27i for pM+G in powers of g^/m^. Other details given in the text. 



Figs. Cni and CSl 

Let us now summarize our findings. We first collect the main results for the optimal approximants (which have 
been presented in detail in Figs. I8ll'2|l and combine them in Figs. 1171 and 1181 thereby focussing on the crucial region 
of low temperatures. For comparison we also include the predictions of lattice calculations taken from Fig. 4(b) of 
Ref. 43] (which includes lattice results of Refs. [23 )• These figures thus present our resummation results when Borel- 
Pade BP[l/2](a(^E)) is apphed to expansion l|SSl) for i?g^", and Fade [0/3](.g|/TOE) to expansion (|77|l forpM+G- The 



central values of the renormalization scales were taken fiE = 2'kT and ^m = Mn 



m 



(0) 



(T) [Eq. 1021)], and the central 



In Fig. El on the other hand, the Ag /mE-term is added separately as the leading order term of expansion (|33|l . i.e., 

z| resummed as P[l/l](a(/im))]- In the 
becomes appreciably weaker, as was 



factorization scale was Ae = (/xeMm)^^^- The EQCD parameters g^ and m| were calculated (resummed) as Fade 
P[l/l]{gs{nm)), at low-energy renormalization scale /xm — ^J-m [^ "^e ^ dsT). In addition, the effect of replacing the 
Pade-resummed P[2/3](a) j3 function by the simple TPS /3 is displayed as well. For the lattice results, we used the 
values of critical temperature T^{nf = 3) « 154 MeV |4J . and Tc{nf = 0) « 267 MeV. The latter value is obtained 
from the result Tc/V^ = 0.629 :24l] and ^ = 425 MeVpl- 

The two Figures IjlTI and lTS)) differ by the different treatment of the Ag -term: In Fig. El the A^ /rriE-term of 
expansion (|27() is added separately as the leading order term of expansion H32() . i.e., as a term proportional to g^ium)- 

as a term proportional to (Se/'^e)'^ [with the EQCD parameters g^ and m 
latter case, the dependence on the low-energy renormalization scale /im = M' 
seen in Fig.El[cf- also the discussion following Eqs. l|^ - l|^ ]. 

In both cases, the lattice results differ from these curves by about 20% and 10% when n/ = 3,0, respectively. 
Taking into account that all lattice results should be taken with an error of 10-15%, we see that our resummed 
results come at low temperatures reasonably close to the lattice results, although our resummations are based only on 
perturbation expansions. Further, in contrast to the various TPS evaluations, the dependence of our results on the 
unknown parameters 6q and (5e and on the high-energy renormalization scale fiE (^ 2t:T) is quite weak as shown in 
Fig. 1171 The dependence on the low-energy renormalization scale /^m — ^J■m (~ '71e(T)) is in Fig. 1171 bv about a factor 
of 2-3 weaker than in the analogous TPS case, and is thus at low temperatures still quite significant. It has its origin 
primarily in the /im-dependence of the QCD term g^{^im) of Eq. 1)321) . In Fig. El the variation of the curves when 
(5g, ^e and /^e are changed, at any given T, are identical to the variations of the corresponding curves of Fig. 1171 i.e., 
quite weak. The dependence on the low-energy renormalization scale /x™ = ^m is in Fig. El appreciably weaker than 
in Fig. El 

Let us recall also, that the specific choices for the approximants which entered in Figs . 1 1 71 and 1 1 81 were based on the 
following reasoning: a TPS at a given order, in principle, allows for various Fade and Borel-Pade approximants. We 
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chose those approximants which give results reasonably stable under the variation of the corresponding renormalization 
scales /^E and /j,m = /imi and of the factorization scale Ag. This physical criterion led us to two possible approximants: 
BP[0/3] and P[0/3] for pm+G: in both cases BP[l/2] for R^'^. They differed from each other significantly in the low 
temperature regime, cf. Fig.|5Jb). In order to further eliminate one of them (BP[0/3]), we had to apply an additional 
physically motivated criterion, namely that the predicted pressure should be smaller than the ideal gas value if the 
temperature T is sufficiently near the critical value Tc- In addition, two numerical arguments were provided which 
also support the elimination of the approxmant BP[0/3] for pm+g- Thereby we ended up with the approximants 
BP[l/2] for i?g^" (^ pe) and P[0/3] for pm+g {=^ Pm+g) which are considered as our "best approximants". 
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FIG. 18: Analogous to Fig. ^ but the A^ V"^E-term of Eq. ED i 



IS now written as 



proportional to {g'^/m^)'^ [Eq. 1331 ] instead of gf 



[Eq. 1321 ]. Changes under the variation of fiM = Mm are included as in Fig. 1171 and (5e = (5g = is kept. The curves vary by the identical 
amounts, at any given T, as in Fig. 1171 under the changes (5e = ±^2 ^'^'^ ^G ~ ^^• 



Since both are of Pade type, we expect that they also have a better convergence behavior than the ordinary 
truncated perturbation theories. This, in fact, can be tested explicitly and the results are presented in Fig. 1191 where 
we demonstrate the behavior of the Pade and Borel-Pade resummed curves, for nf = 3, when the order of expansions 
(|38|l and H27|l for R'^^ and pm+g increases, and compare them with the corresponding TPS evaluations. For all the 
curves here, including the TPS curves, the EQCD parameters m| and g\ were evaluated as P[l/l](a(/im)), as in 
Figs. m and UTUHl The 0{g^) P-curve uses P[l/l](a(^E)) for i?|f" and P[0/l](.g|/mE) for pm+g of expansion |23; 
the 0{g^) P-curve uses P[l/1] for i?f " and P[0/2] for pm+g- The 0{g^) (P+BP)-curves use BP[l/2] for i?f " and 
P[0/3] for Pm+g ~ the lower curve is the central curve of Fig. El the upper is the central curve of Fig. ^J In Fig. 1191 
we see that the TPS results change strongly when the order of the TPS is increased; on the contrary the resummed 
results suffer weaker changes, although a clear convergence at T < 1 GeV still cannot be seen at these orders. 

The crucial point of our approach was to treat separately the short-distance (pe) and long-distance (?Tig, ^g, pm+g) 
quantities - using in them the renormalization (and factorization) scales which correspond roughly to the physical 
scales of the considered quantities, and then performing either Pade or Borel-Pade resummation of each quantity. 
In this way we arrive at predictions for pressure which have weak dependence on the unknown parameters Je and 
(5g, and on the high-energy renormalization scale [i^ {^ 2ttT) of pb- On the other hand, the dependence of the 
predicted pressure on the renormalization scale fim = fJ-M (~ tie) of Pm+g is still large at T < 1 GeV when the 
EQCD Ag /wE-term in the EQCD expansion (|27|l is added separately as a power of the QCD coupling parameter 

gsilJ-m) [Eq. 1321) ]. This /x„i-dependence becomes significantly weaker when the EQCD Ag /?TiE-term is instead added 
separately as a power of the first EQCD coupling parameter g|/m.E [Eq. I|33|l]. Further, the resummed results change 
reasonably slowly when the order of the TPS's, on which they are based, is increased. 

On the other hand, simple TPS evaluations do not yield reasonable results for the long-distance quantities at low 
temperatures, because the TPS's at such low energy scales have strong dependence on the unknown parameters and 
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FIG. 19: The convergence behavior of our (Borel-)Pade curves for p/pidcal ^s a function of temperature, at rif = 3, for orders 0(gf ), 
0{gg) and 0(51^), and the behavior of the coresponding TPS's. The upper P+BP ©(gf ) curve is for the case when the A^ /m^-term of 
Eq. 1271 is now written as proportional to (gg/rriE)'' instead of gf . Details are given in the text. 

on the renormalization scales, and change strongly when the order of the TPS's is increased. 
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APPENDIX A: PADE AND BOREL-PADE APPROXIMANTS 

In this Appendix, we describe the procedure we used to obtain the Fade and Borel-Pade approximants on the bases 
of perturbation expansions (PS|) . (|25J|, l|M|) . (P7|) and (j^Hl- In expansions (P5|l . if^ and (|55|) . the expansion parameter 
was a(/l) = [gs (/l)/27r] . The canonical expansions of this type have the form: 

Sia) = all + Y,rnaA . (Al) 

Fade approximant of order {N,M) is P[N/M]s{a) = Pn {a) / Pm {a) , i.e., the ratio of polynomials PAr(a) and Pm{o) 
of order N and M , respectively. The coefficients of these polynomials are determined by the requirement that the 



re-expansion of the approximant in powers of a reproduce the terms up to (and including) 
(jAlj) . For example, 

F[l/l](a) = 



-.N+M 



(1 - r^a) 



in the expansion 



(A2) 



The Borcl-Pade approximation BP[A^/M]s(a) is constructed from expansion IjAip by applying Fade V[N /M]Bs{b) to 
the expansion of the Borcl-transform 






(A3) 
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and then performing the Borel integration 

BP[N/M]s{a) ^^ r dh exp f--^) P[iV/Af]B,(6) . (A4) 

Po Jo V Poa/ 

In principle, the integration here is along the positive 6-axis. In order to avoid any possible poles on the positive 
&-axis, the integration is in general performed along a ray in the b plane, say: b — rexp(— i0), with </> ^ arbitrary 
and fixed; then r = |6| is integrated from zero to infinity, and the real part of the result is taken - cf. Ref. jsj- This 
leads, by Cauchy theorem, to the Principal Value (PV) prescription for integral (|A4I) . 
On the other hand, expansions (|5H) and l|77|) for pm+g have the following form: 

oo 
n=l 

where g = gsifJ-u) or g = g|/mE- Pade approximant P[N/M]T{g) = PN{g)/PM{g) was in this case simply the Pade 
applied to expansion (|A5p . The Borel-Pade approximation BP[A^/M]t(5) was constructed from expansion (|A5|I by 
applying Pade V[N/M]~ [z] to the (expansion) of the following Borel transform of T: 

BT(z) = l + f;^z", (A6) 

n—l 

and then performing the corresponding Borel integration 

BP[N/M]T{g) ^- I dz exp (--\ V[N/M]~^{z) . (A7) 

Again, in order to avoid the possible poles at z > 0, the integration is performed along a ray in the z-plane and the 
real part is taken. 
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